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In a recent paper [M. A. J. H. Ossendrijver, A&A 359, 364 (2000)] numerical simulations of a 2D 
mean-field model where shown to produce grand minima, typical of the long-term behavior of solar 
magnetic activity. The model consisted of dynamo that features an a effect based on the buoyancy 
instability of magnetic fiuxtubes, which gives rise to the switching back and forth from grand minima 
to "regular" solar behavior. In this Letter, we report evidence from a time-series analysis of the 
model for the presence of crisis-induced intermittency due to attractor-widening. We support this 
finding by showing that the average duration of the minima, (r) , follows the theoretically predicted 
scaling (r) ~ {Csa — C|^) ', where Csa is the bifurcation parameter of interest, together with other 
statistical evidence. As far as we are aware, this is the first time concrete and detailed evidence 
has been produced for the occurrence of this type of crisis-induced intermittency - due to attractor 
widening - for such dynamo models. 



The records of past solar magnetic activity reveal an 
outstanding phenomenon referred to as grand minima. 
During the latest of such intervals, the so-called Maun- 
der minimum (1645-1715), sunspots were virtually ab- 
sent. This and earlier grand minima are clearly visible in 
the records of cosmogenic isotopes, e.g. ^^C and ^"Be [g|). 
Timing and duration of the known grand minima are ir- 
regular. Solar variability is also apparent e.g. in length 
and amplitude variations of the 11-year sunspot cycle. 
The spectral and statistical properties of solar variability 
are still not well-known. Although there is some evidence 
for various modulations of the solar cycle H , it remains 
unclear, due to the lack of a sufficiently long and accu- 
rate data set, whether they are truly periodic rather than 
chaotic or even purely stochastic Q] . 

Intermittency has been proposed p] as one of the pos- 
sible scenarios for the underlying mechanism behind the 
occurrence of grand minima. Several dynamo models 
have been show to exhibit quasiperiodic or chaotic in- 
termittent behavior 0. In |jl|] a case was made for a 
stochastically driven 2D mean-field dynamo model show- 
ing grand minima [[7|. Earlier, the authors in g produced 
grand minima using a ID model, in which the flux injec- 
tions were treated as an additive random source term. 
The main advantage of a 2D model is that the radial 
structure is resolved, so that the flux injections can be 
treated more realistically through the advection term. Of 
course, mean-field theory cannot replace full MHD cal- 
culations ||] and is at best capable of capturing the most 
important aspects of solar dynamo action. The purpose 
of the calculations was to illustrate that grand minima 
are an inherent feature of a solar dynamo based on mag- 
netic flux tubes. 

In this Letter, we show concrete evidence that the 
switching back and forth from grand minima to the "reg- 
ular" 22-year cycle is a manifestation of a known dy- 
namical process called attractor-widening, resulting in 
crisis-induced intermittency fin] . This type of intermit- 
tency has been found in several experimental and nu- 



merical studies pTl. Other types of crisis pOJ have also 
been concretely demonstrated to exist in dynamo models 



I. DYNAMO MODEL 

The induction equation for the mean magnetic field B 
in the first- order smoothing approximation is given by 



n T3 

-^- = V x{uxB + aB-(iV xB]. (1) 

at 



All mean quantities are defined here as longitudinal av- 
erages, so that the mean field is the axisymmetric com- 
ponent of the actual field. The a term parametrises the 
effect of helical motions. In the solar dynamo, its main 
effect is to create a poloidal field from a toroidal field. 
The total a effect consists of a buoyancy-driven term 
in the overshoot layer, which exists only if \B\ exceeds 
a threshold value, and a small-scale kinematic a effect 
in the convection zone, which is modeled as a random 
forcing term. The flow u consists of differential rota- 
tion (chosen to model helioseismological inversion results 
Pjl), one updraft at a fixed low latitude, and downdrafts, 
with a finite lifetime, placed at random chosen latitudes. 

For a detailed description of the model, the reader is re- 
ferred to [|l| ('model B'). Its main features are as follows. 
The dynamo equation (n^) is solved in a spherical shell 
between r = O.5i?0 and r = i?Q, consisting of a layer 
with convective overshooting between ri = O.6i?0 and 
r2 — O.7i?0, and a convection zone above it. The mag- 
netic field strength is expressed in units of i?o, the thresh- 
old value for the buoyancy instability. Time is measured 
in terms of i?^//3o, the magnetic diffusion time for the 
convection zone (see below). 



II. RESULTS 

The numerical integration is carried out using an im- 
plicit 2D code developed by D. Schmitt and T. Prautzsch 
of Gottingen. The grid size is set to 61 points in the ra- 
dial direction, and 51 in the latitudinal direction. 

Fig. |l| shows a butterfly diagram for the mean toroidal 
magnetic field of a solution with grand minima (model B, 
Fig. 4 in [Q). The numerical solution features intervals 
of ordinary cyclic dynamo action interrupted by intervals 
of irregular duration without cycles, reminiscent of the 
grand minima of solar activity. The ordinary solar cycle 
is characterized by belts of magnetic activity on both 
hemispheres, which have opposite polarity, and migrate 
towards the equator during the course of a cycle. At the 
solar minimum, the old belt disappears, and a new belt of 
opposite polarity is formed at higher latitudes. Note that 
the present schematic model is not designed to reproduce 
all known detailed features of the solar cycle but some 
of the most important aspects, particularly the regular 
butterfly diagram and the existence of grand minima. 




FIG. 1: Butterfly diagram of the toroidal field ai r — 
O.6257?0 for parameters C^ = 2-10'', Cd = -0.05, Cc2 ^ -1, 
Csc = 100, Cd = -0.1, and Cb = 0.014 (from Ossendrijver § 
- model B) . Cu measures the strength of the differential rota- 
tion, Cai the strength of the a effect on the overshoot layer, 
Csa the strength of the mean a effect in the convection zone, 
Ca2 the strength of the fluctuating a effect in the convection 
zone, Cd the strength of the downdrafts and Ct the strength 
of the updrafts. Time is in years. 

Fig. H shows the evolution of the solutions as one pa- 
rameter of interest, Csa, is changed. The shutdown peri- 
ods correpond to grand minima and the active period cor- 
rcpond to the (noisy) solar cycle. As can be easily seen, 
if C'sa is raised, the bursting becomes more and more 
frequent. The critical parameter value Cg^ for which the 
smaller attractor, corresponding to the minima, starts to 
burst was found to be around C^ ^ 59. This is indicative 
of the existence of some form of crisis, which we substan- 
tiate below. According to the theory |jlfl, around the 
critical parameter value C^ ^ 59 the attractor collides 
in phase space with the stable manifold of some unstable 
periodic orbit. The excursions or burst phases seen in 
Fig. result in the widening of the attractor after this 
crisis. 

Phenomenologically, the time series could be inter- 
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FIG. 2: Times series of the toroidal field, in units of Bo, at 
r — 0.625Rq and 6 — 0.477r. Except the bifurcation param- 
eter Csa, all parameters are as in Fig. hi As Cfo is increased 
beyond the critical value, C^ ~ 59, the frequency of the shut- 
downs (grand minima) decreases rapidly. 



preted in terms of various different types of intermit- 
tency. In particular, the time series could be mislead as 
solutions of systems exhibiting on-off intermittency [15|, 
or Type I/II/III Pomeau-Manneville intermittency [16|. 
However, we have eliminated this possibility by analysing 
both the power spectra (see Fig. y) and the probability 
distribution of the shutdown phases (see Fig. Q) . 

It has been shown [|l^ that the corresponding power 
spectra for Type I intermitency have a broad-band fea- 
ture whose shape obeys approximately the inverse-power 
law, P{f) ^ /~^, for / > fg, where /^ is the saturation 
frequency. Our power spectra do not seem to show this 
property. Neither do they conform to the theoretical ex- 
pectation that the power spectra for an on-off process fol- 
lows a law P{f) ~ f~^^^ II3I over an intermediate range 
of frequencies. Furthermore the probability distribution 
of shutdown phases (Fig. Q) does not follow the scaling 
P{n) ~ n~2j as predicted for on-off intermittency |l^ 
but the theoretically predicted scahng ijl^l 



r/{T 



(2) 



Finally we confirm the agreement with the theoreti- 
cal predictions for "crisis-induced intermittency" by cal- 
culating the average time (r) between the bursts which 
obeys the scaling law |1^ 



{t) r^ ICSa - QJ 



(3) 



which we confirmed for our model, as can be seen in 
Fig. |. The value of 7 is found to be 7 = 2.19 ± 0.04. 
This further reinforces the conclusion that the type of 
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FIG. 3: Power spectra of the time series in Fig. m with 
Cfa ~ 100. Note the comparison with known scalings of power 
spectra for solutions of systems showing Type I intermittency 
or on-off intermittency. It shows that the time series cannot, 
in principle, be due to these types of intermittency. 



10" 



10" 



Pi 



10" 



10" 




0.1 



1.0 



10.0 



FIG. 4: Scaling of the probability distribution of the du- 
ration of shutdown phases for the time series in Fig. Q, with 
C<L ~ 100- Good agreement is found for both the toroidal and 
poloidal (not shown) numerical obtained distribution with the 



known theoretical scaling, r ~ e 
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intermittency observed here is neither Type I (for which 
7 = — ^), Type II or III (for which 7 = —1) or on-off 
intermittency (for which 7 = — 1). 

The 7 coefficient can be calculated also analytically, as 
shown in ||l^ . The method involves calculating the stable 
and unstable manifolds of the unstable orbit mediating 
the crisis. However, because our attractor seems to be of 
a high dimension, the method of using the time series to 
reconstruct the attractor and analyse Poincare sections 
of such embeddings was not useful in identifying the me- 



diating unstable orbit. The value of 7 seems to confirm 
that our attractor is indeed of a high dimension, given the 
theoretically expected range {D — l)/2 < 7 < (-D + l)/2 
[jTo). This would imply D = 4. 
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FIG. 5: Scaling of the average times of shutdown (minima) 
(r) as a function of {Csa — Csa) for crisis-induced intermit- 
tency due to attractor widening for model {}^. The slope is 
found to be 7 = 2.19 ± 0.04 and the error bars are ±(r)/y^, 
where n is the number of minima found for each parameter 
value. 



III. DISCUSSION 

We have obtained concrete evidence, in terms of time 
series signatures, power spectra and dynamical scalings, 
to demonstrate concretely the presence of crisis-induced 
intermittency due to attractor widening in a mean-field 
PDE dynamo model. Despite the presence of simplifi- 
cations in these models, this is of potential importance 
since it shows the occurrence of another type of intermit- 
tency (in addition to Type I Pomeau-Manneville, attrac- 
tor merging intermittency [O and in-out intermittency 
|l9| recently discovered) in these models. This may be 
taken as an possible indication that more than one type 
of intermittency may occur in solar and stellar dynamos 
pO| . Given that the underlying model showed several 
important features of solar activity, such as minima of 
activity and a realistic butterfiy diagram, this suggest 
that intermittency, and in particular the form showed 
here, may be the underlying nonlinear dynamical pro- 
cess responsible for the observed solar magnetic activity. 
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